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THEIR GEOMETRIC AND RADIOMETRIC ACCURACIES 


by 
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ABSTRACT 

Some Landsat Thematic Mapper data collected over the eastern United States 
were analyzed for their whole-scene geometric accuracy, band-to-band 
registration and radiometric accuracy. Band-ratio images were created for a 
part of one scene in order to assess the capability of mapping geologic units 
with contrasting spectral properties. Systematic errors were found in the 
geometric accuracy of whole scenes, part of which were attributable to the 
film-writing device used to record the images to film. Band-to-band 
registration showed that bands 1 through 4 were registered to within one 
pixel. Likewise, bands 5 and 7 also were registered to within one pixel. 
However, bands 5 and 7 were misregistered with bands 1 through 4 by 1 to 2 
pixels. 3and 6 was misregistered by 4 pixels to bands 1 through 4. 

Radiometric analysis indicated two kinds of banding, a modulo-16 stripping and 
an alternate light-dark group of 16 scanlines. A color-ratio composite image 
consisting of TM band ratios 3/4, 5/2, and 5/7 showed limonitic clay-rich 
soils, limonitic clay-poor soils, and nonlimonitic materials as distinctly 
different colors on the image. 
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A PRELIMINARY EVALUATION OF LANDSAT-4 THEMATIC MAPPER DATA 

FOR 

THEIR GEOMETRIC AND RADIOMETRIC ACCURACIES 


INTRODUCTION 


This report describes results of some preliminary analyses of Landsat-4 
Thematic Mapper data for the NASA Landsat Image Quality Analysis program 
tLIDQUA). The work is being done under interagency agreement S-10407-C 
between the U.S. Geological Survey and NASA-Goddard Space Flight Center vGSFC). 

Landsat-4 TM scenes for Washington, D.C. (40109-15140, November 0, 1930' 

Macon, Georgia (.40050-15333, September 4, 1930) and Cape Canaveral, Florida 
v4003o-l4014 , August 00, 1930) have been examined to determine their geometric 
and radiometric accuracy. In addition, part3 of these scenes are al30 being 
analyzed to determine the ability to identify specific rock types with the 
added near-infrared TM bands at 1.6 and 0.0 urn. 

GEOMETRIC ANALYSIS 


In crier to assess mole scene geometric accuracy, each of the seven TM b\r.ds 
for the geometrically corrected vGSFC ?-tape digital format) Washington, 0.'. 
scene were contrast stretched ano recorded on film using an Optronics P-lcCO 
film writer. Thirty-eight well-defined control points (.figs. la and cl were 
visually identified by comparing 1 : 0 4 .CCO-scale topograpnio taps with TM oa.nd 
Tender and 'll "re* stated mat th-» * \r.bir.eo affect cf mar-pcsiticr. 

>rror (.for maps meeting National Map Standards' ana control-point 
identification error may yield a Root 'car. Square Error (RMSE of 10 r . a 
point marking device was usee to drill d0 micrometer circular aoles at me 
control points in me emulsion of the band 3 image. The control points “.men 


.ere .ere trar.s: erreo am 


- * k. 


' onto .mages of bands 1. 5. and 


oinocuiar vision and tme tnllec holes on card 3 as the '’master". Therefore, 
ar.\ errors arising luring the transfer of me control points from the mac : c 
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the hand 3 image also are present in the other bands. Cartesian coordinates 
of each control point on each band then were measured by means of a 
photogrammetric comparator whose resolution is one micrometer. 

Universal Transverse Mercator (UTM) projection coordinates for each of the 
control points were measured from the topographic maps and recorded to the 
nearest 5 m using a metric coordinate reader (a transparent plastic 
template) . The UTM coordinates were transformed numerically to 
latitude-longitude and then to Space Oblique Mercator (SOM) coordinates. 

The coefficients of two transformations, the similarity and the affine, were 
calculated by the least-squares method by using the measured comparator 
coordinates, and their respective UTM and SOM coordinates. The similarity 
transformation can be interpreted as a rotation, single change in scale and 
two translations: 


x' : ax + by + x 0 a = X cosO 

y ' = -bx ♦ ay + y Q b = X sinO 

when X is the scale change, 0 is the rotation, and x 0 , y 0 are the 
translations . 

The affine transformation is a six-parameter transformation permitting two 
mutually perpendicular scale changes, two rotations (one defining the 
direction of the scale changes), and two translations: 

x* = ax by + x 0 

y' = cx dy + y 0 

Table 1 shows the vector residuals resulting from each transformation for the 
four TM bands measured to date. 


Table 1 

ROOT MEAN SQUARE ERROR ( RMSE ) OF VECTOR RESIDUALS 
FROM CONTROL MEASUREMENT 


UTM ground coordinates SOM ground coordinates 


TM band 

Similarity 

Affine 

Similarity 

Affine 

1 

72 a 

45 m 

70 m 

4 3 3 

3 

70 m 

45 m 

67 m 

42 m 

5 

63 a 

43 a 

0 0 -w 

4q 71 

7 

69 m 

45 2 

65 m 


In all cases , 

the SOM projection fits 

the data sli 

■fthtly better, i 

.e. r.as a 

smaller RMSE, 

tr.an one UTM projection, 

out little 

^ ^ a r» a n £ o ^ ^ ^ 

hate exists 


between a scene recorded on film in eit.tsr of tr.e tv 
difference prooaoly will be smallest for tncse areas 
coincides with a UTM central meridian; at these loca 


.n which tr.e scene cer.te 
cas the nigher order 


terms of the two projections are minimal. 
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Table 1 shews that the application of the affine transformation significantly 
reduces the root mean square error in all cases. One possible cause may be 
the film-writing device. A highly significant sxewness had been observed in 
1930 from tests of grids printed from the same device as used in this study 
(N. Falcone, U.S. Geological Survey, personal communication). At that time, 
it was hypothesized that it possibly was caused by lack of parallelism between 
the axis of the rotary drum and tne axis of travel of the light source. In 
order to verify the I960 findings, a new six by six grid, having 4 cm spacing, 
was written on film using a 50 micrometer aperture. The RMSE values are given 
in Table 2 along with those of the 1980 measurements: 


Table 2 

ROOT MEAN SQUARE ERROR OF VECTOR RESIDUALS 
FROM GRID MEASUREMENTS 

Similarity Affine 

1980 Grid 136 micrometers 75 micrometers 

1983 Grid 117 micrometers 21 micrometers 

Although the values for the 1983 grid are smaller, a significant affine error 
is still indicated. Vector plots of residuals from the similarity 
transformation display an extremely high degree of correlation between the two 
sets of grids. Figures 2 and 3 are vector plots of residuals for the 1983 
grids. Figure 2 indicates a general elongation inclined at 45° to the 
direction of the lead screw which moves the light source and a compression at 
90° to the prior direction. Figure 3 is a plot of residuals remaining after 
the removal of first order errors. One can surmise, from examination of the 
residuals, that at least a third order polynomial is necessary to reduce the 
residuals to the same order of magnitude as the measuring errors 
(approximately 3 micrometers) . 

In order to assess the difference in shape between the UTM and SGM map 
projections, their respective sets of coordinates of the control points were 
fit together by the least-squares method. The similarity transformation 
results in a RMSE of 11 . 3 m for the 38 control points, whereas the affine 
transformation results in a value of 11.4 m. Because both projections are 
conformal, it is not surprising that the affine transformation does not reduce 
the residuals. 


Sand-to-band registration was checked for the Washington, D.C. scene viewing 
several small parts of the scene on an interactive video display device. 
Several strongly contrasting targets, such as airport runways and small pcr.ds, 
were used 00 visually determine the offsets. By rapidly toggling between the 
TM bands on the interactive system it was found that bands 1 through 4 were 
coregisterec nearly perfectly (fig. -a). Using channels 1 through 4 35 
references, bands 5 and 7 were found to be displaced 1 to 2 pixels westward 
(fig. 4b). Band 0 was displaced 4 pixels eastward (fig. 4c), wnich 
effectively means it is displaced by one of its 120 n resolution units. 



RADIOMETRIC ANALYSIS 


^piQ^AL PAG'- : -■ 

OF POOR QUALITY 

A geometrically uncorrected version (GSFC A-tape digital format) of the Macon, 
Georgia scene was examined to determine some of the radiometric 
characteristics of the data. This included analysis of the 16 individual 
detectors in bands 1 through 5 and 7 and each of 4 detectors of channel 6. 
Systematic gaps were found within the limited range of digital numbers 
recorded for the scene. These can be attributed to the integer truncation of 
decimal numbers involved in the calibration procedure during production of the 
computer compatible tapes. 

The potential of sensor striping caused by the separate gains and offsets 
required for each of the 16 sensors of the visible and near-infrared TM bands 
was investigated using the A-tape of the Macon scene. The four sensor thermal 
infrared band also was examined. Means and standard deviations were 
calculated on a modulo-16 basis. The results are summarized in table 3* 


Table 3 

COMPARISON OF DIGITAL NUMBER STATISTICS 
FOR MODULO-16 STRIPING OF TM DATA 


TM Channel 

1 

2 

3 

4_ 

5 

6 

7 

range of 

85.5 

34.2 

30.6 

74.5 

66.8 

126.0 

22.7 

sensor means 

86.4 

34.8 

31.2 

74.9 

67.1 

126.7 

23.1 

difference 

0.9 

0.6 

0.5 

0.4 

0.3 

0.7 

0.4 

range of 
standard 

17.649 

9.067 

11.295 

16 .665 

23.449 

15.527 

11.679 

deviations 

18.076 

9.504 

11.661 

17.026 

23-931 

15.980 

12.271 

difference 

0.427 

0.437 

0.366 

0.361 

0.482 

0.453 

0.592 


The range of means is typically less than that of Landsat MSS values when the 
MSS data are scaled into 8-bit space. Enhanced-contrast images of parts of 
the Macon scene were recorded onto film and visually examined for evidence of 
striping. Striping is not evident in channels having sensors with a range of 
mean differences less than 0.4 between their high and low values. Striping 
becomes apparent when difference values reach 0.5 and becomes objectionable 
when values are greater than G.6. 
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of 28. 5 m. The relationship should hold for at least a snail contiguous group 
of scanlines. Two areas along the sane groups of 1~ scanlines were sanpied 
over relatively homogeneous bodies of water. Statistics were contrasted for 
waters of Lake Weir, Florida, on the western side of the image, against the 
Atlantic Ocean waters located off Edgewater, Florida, near the east edge of 
the image (Tables 4 and 5). 


TABLE 4 

COMPARISON OF DIGITAL NUMBER STATISTICS, TWO HOMOGENEOUS 
TER EODIES ALONG THE SAME, BUT OPPOSITE ENDS OF GROUPS OF 

IT SCANLINES 


Lake water 


Scanlines 

TM Band 

1 

2 

3 

4 

5 

6 

7 

2809-0825 

Mean 


TS.97 

29.07 

23-01 

17.22 

9.41 

131.51 

4.24 


Std. 

Dev 

. 1.44 

0.63 

0.66 

0.47 

0.63 

0.50 

1.08 

2826-2842 

Mean 


T9. 92 

29.33 

23.34 

17.36 

9.41 

131.63 

4.27 


Std. 

Dev 

. 1.46 

0.61 

0.73 

0.49 

0.73 

0.47 

0.96 

2343-2859 

Mean 


T3.T4 

29.00 

23.00 

17.15 

9.36 

131.63 

4.21 


Std. 

Dev 

. 1.42 

0.68 

0.70 

0.43 

0.79 

0.47 

1.05 

2860-28T6 

Mean 


79.24 

29.13 

23.22 

17.37 

9.40 

131.72 

4.28 


Std. 

Dev 

. 1.54 

0.68 

0.74 

0.52 

0.79 

1.45 

1.04 

28TT-2893 

Mean 


73.57 

28.32 

22.99 

17.07 

9.38 

131.71 

4.31 


Std. 

Dev 

. 1.60 

0.69 

0.72 

0.45 

.74 

0.46 

1.04 




Sample 

size for 

each group 

of scanlines 

= 272 




TABLE 5 


COMPARISON OF DIGITAL NUMBER STATISTICS, TWO HOMOGENEOUS 
WATER BODIES ALONG THE SAME , BUT OPPOSITE ENDS OF GROUPS OF 

IT SCANLINES 


Ocean water 
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5otn these areas exhibited banding in tr.e image data. The beginning of the 
groups of 17 scanlines was determined by using a harshly contrast-enhanced TM 
band 3 image and picking the oeginning scanline on an interactive video 
analysis system. Examination of table 5 shows that the ocean waters have a 
consistent variation in mean values for groups of 17 adjacent scanlines, 
beginning with all channels naving high mean values for the first group of 17 
lines (2809-2325) , and proceeding in an alternating low-high fashion for 
following groups of 17 lines. Digital values for the laxe waters (table 4) 
show an opposite trend. The same initial group of scanlines shows 
consistently low values and proceeds in an alternating high-low fashion for 
following groups of 17 scanlines. Differences between adjacent groups of 
scanlines are greatest in TM channels 1 through 4, and appear negligible in TM 
band 5-7. The change in. relative mean values across the image for adjacent 
groups of scanlines suggests an inherent and systematic change occurring, 
perhaps due to some decay factor imposed on the gain or offset values for 
opposite sweeps of the pushbroom scanner. Some crossover point must exist in 
the data, presumably near the center of the scanline, wnere digital numbers 
for opposite sweeps would be of equal value. It also suggests that this type 
of banding could not be corrected by a modulo-16 or -32 histogram 
equalization, even in the geometrically uncorrected A-tape data, unless 
correction factors which varied as a function of position along a scanline, 
were applied along groups of 32 scanlines. 

PRELIMINARY SPECTRAL OBSERVATIONS 


Specific band ratios were used to determine their effectiveness in defining 
vegetation and rock mineralogy for the Washington, D.C. scene. The TM 3 /4 
band ratio was selected to depict vegetation. Vegetation appears black to 
dark gray in this ratio because of its relatively low reflectance in band 3 
due to the presence of chlorophyl absorption bands, and its high reflectance 
in band 4 due to a lack of chlorophyl aDsorption (Knipling, 1970). Geologic 
materials lack the strong contrast between TM bands 3 and 4 and thus would 
appear light gray to white in the ratio. The 5/2 band ratio was used to 
depict limonitic rocks. Limonitic rocks have low reflectances in bands 1 and 
2 compared with band 3 due to a strong ferric-iron absorption band centered in 
the ultraviolet, but whose effects reach into the short wavelength visible 
part of the spectrum (Hunt and others, 1971; Hunt and Ashley, 1979). 
Nonlimonitic rocks and soils lack the absorption band in the ultraviolet, 
hence their spectral curves do not show the sharp drop between band 3 and 
bands 1 and 2. Most geologic materials have no absorption bands in TM band 5 
and are quite bright in this region of the spectrum. Hence, a 5/2 band ratio 


image will appear relatively bright for limonitic rocks, whereas nonlimonitic 


!oks will appear dark. 


Rocks which have minerals that contain hydroxyl or carbonate anions in 
crystal lattices, have absorption bands v Hur.t and Salisbury, 15~1; Hunt 
Ashley, 1579 in the bandpass of TM oar.a ~. Hydroxyl-bearing minerals, 
containing either Al-C-H or Mg-C-H bands, include: 1 clay minerals, 2 

sheet silicates suer, as the serpentine minerals, some of the micas, and 
among others, anc 3) sulfate minerals such as adr.ite and jarosite, vhi 
often are found as products of hydrothermal alteration “recesses . Tart 
minerals induce calcite anc dolomite, wrier, are common constituents m 
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sedimentary rocks limestone and dolomite. Rocks containing these anions will 
have relatively low reflectance values in TM hand 7, whereas rocks lacking 
these anions will have relatively higher reflectances. A 5/7 hand ratio was 
evaluated to test the effectiveness of the 1.6- and 2.2-um bands in detecting 
rocks containing hydroxyl- or carbonate-bearing ainerals. Areas containing 
rocks which have absorption bands in the 2.2-um region will appear lignt gray 
to white in the 5/7 ratio, whereas rocks lacking absorption bands in this 
region will appear dark gray to black. Similar bands were selected in prior 
studies to map clay-rich rocks using airborne Thematic Mapper Simulator (TM S) 
scanners (Podwysocki and Segal, in press; Rowan and Kahle, 1932; Ashley and 
Abrams, 1980). Preliminary results show that clay minerals containing 
absorption bands at 2.2 um, which are centered in TM band 7, can be readily 
detected when the band ratio images are suitably contrast-enhanced. 

Figure 5a is a false-color infrared image which shows a 256 by 256 pixel area 
from the Washington, D.C. scene. A road construction project has laid bare 
clay-rich soils containing ferric iron oxide minerals. A color ratio 
composite image (fig. 5b) shows the clear barren area as white, which in the 
color assignments made to each ratio, indicates both ferric iron oxide and 
clay minerals. Ferric-iron-rich areas lacking clay minerals are portrayed as 
a cyan color. Vegetation appears red to magenta, depending upon its type and 
vigor. Spectrally flat areas (lacking both ferric iron and clay absorption 
bands) are green to yellow-green. 

SUMMARY 


Lar.dsat-4 Thematic Mapper digital data were examined for several scenes in the 
eastern United States to characterize their geometric and radiometric 
accuracies. A small portion of one of the scenes also was processed using 
band ratioing in order to test the radiometric characteristics for mapping of 
spectrally dissimilar geologic materials. 


Whole-scene accuracy was assessed by comparison of 38 image and topographic 
control points for four bands of a TM image of the Washington, D.C. area. The 
images were compared with both Universal Transverse Mercator and Space Oblique 
Mercator projections. Both similarity and affine transformations were applied 
using the least-squares method in order to fit the image control points to the 
map coordinates. Results show that the affine transformation provides a 
better fit of the image to the map projections, decreasing the root mean 
square error from an average of about 70 m for the similarity transformation 
to about 45 m for the affine transformation for both projections. Results 
show that little difference exists between the two projections for this TM 
scene. .Analysis of the geometry of the film-writing device that was -sec to 
produce tre TM images suggests a systematic error in the device, causing an 
elongation along one diagonal of the image and a foreshortening along tr.e 
oppcs-te diagonal. This distortion in the film-writer device accounts for 
most residual error found in tr.e similarity transformation of the TM image 
data. Hence, geometric errors in the TM data itself, are of some value lass 
than or equal to tr.e RMSE values derived from tr.e affine transformation. 


3ar.d-to-oar.d registration for portio 
subpixel misregistration in TM oar.os 
displaced between 1 to 2 pixels from 
througn 4 and oar.d : is displaced by 



their corresponding points in the bands 1 
3 to J pixels from bands 1 through 4 . 
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Analysis cf the radiometric quality shows at least two types of scanline- 
relatec banding. Radiometric stripping occurs on a mcdulo-16 basis and is 
relatea to small shifts in gains and/or bias between each of the It detectors 
in each band. Stripping is most noticeable in band 1 and oecomes less 
noticeable with increasing wavelength in the visible and near-infrared bands. 
The thermal band six also has relatively pronounced banding. Another type of 
banding appears as alternating groups of 16 lighter and darker scanlines. The 
effect is noticed over homogenous areas such as water bodies and is most 
pronounced along the edges of the scenes. The effect changes along the length 
of the groups of 16 scanlines. Along one edge of the scene, average values of 
the 16 scanlines alternate in a high-low fashion, whereas on the opposite 3ide 
of the scene, values alternate in a low-high fashion for the same scanlines. 
The severity of tne banding is worse in band 1 and becomes less and less 
noticeable towarcs the near infrared. TM band 6 is not affected by the 
phenomenon. 

TM bands were ratioed in order to test the radiometric quality of the data for 
mapping geologic units with contrasting spectral characteristics . A 3/4 band 
ratio was processed to distinguish between vegetation and other materials. A 
5/2 ratio was used to distinguish bwtween limonitic and nonlimonitic rocks and 
a 5/7 ratio was used to detect geologic materials which contained minerals 
having the hydroxyl or carbonate anion. Results showed that clay and 
limonite-rich rocks could be distinguished from limonitic rocks lacking 
clays. Rocks lacking both limonite and clays could also be distinguished. 
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Figure lb - Schematic map showing the distribution of 38 control points chosen 
for this_ study . 
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DIRECTION OF 1 EAD SCREW 
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A TM false-color infrared composite of bands U, 3 and 2 
projected as red, green, and blue, respectively. The 
bright white linear area in the northeast quadrant (C) is 
an area of disturbed ground along a highway right-of-way 
containing clay-rich limonitic soils. Distance across the 
bottom edge of the image is approximately 7*3 km. 
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Figures 4 


Color-composite images (64 by 5u 
jectec or. ire CRT ns rea , green, 
for a portion of Dulles Internal! 
anc sma.j. * a e were use*! to study 
tion. Distance along bottom edge 
approximately 1.3 xm. 


pixels) of TM bancs pro- 
ar.d blue, respectively, 
oral Airport. The runways 
band-to-band registra- 
of the image is 



Figure 4a - A composite image of bands 4, 3, and 2 . Misregister 
appears to be at the subpixei level. Note that the 
brightest portion of the airport parking ramp (A) is 
centered in the middle of the ramp (cf. fig. 4b). 



A TM false-color infrared composite of bands 4, 3 and 
projected as red, green, and blue, respectively. The 
bright white linear area in the northeast quadrant (C 
an area of disturbed ground along a highway right-of- 1 
containing clay-rich limonitic soils. Distance acres 
bottom edge of the image is approximately 7.3 km. 




